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Summary. In this article I will review some basic results on elliptic boundary value 
problems with applications to General Relativity. 



1 Introduction 



Elliptic problems appear naturally in physics mainly in two situations: as 
00 . equations which describe equilibrium (for example, stationary solutions in 

■ General Relativity) and as constraints for the evolutions equations (for ex- 

ample, constraint equations in Electromagnetism and General Relativity). In 
addition, in General Relativity they appear often as gauge conditions for the 
evolutions equations. 

""^ ■ The model for all elliptic equations is the Laplace equation. Let us consider 

the Dirichlet boundary value problem for this equation 

*V Au = f on Q, u = g on dQ, (1) 

where fl is a bounded, smooth, domain in M™ with boundary df2] /, g are 
smooth functions and A is the Laplace operator in M™. 

It is a well known result that for every source / and every boundary 
' value g there exist a unique, smooth, solution u of JQ). We would like to 

generalize equations Q for more general operators and more general boundary 
conditions. 

The first step in this generalization is given by the Neumann problem 

Au = f on ]?, n L d i u = on <9i?, (2) 

where n l is the outward unit normal to df2, the index i takes values i — 
1, • ■ • , n and di denotes partial derivative with respect to the l n coordinate 

There exist two main differences between the Neumann and the Dirichlet 
problem: (i) The solution to the Neumann problem is not unique, for a given 
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solution we can add a constant and obtain a new solution. Moreover, the 
constants are the only solutions of the homogeneous problem 

Au = on f2, n l diu = on dtt. (3) 

To see this, we multiply @ by u and use the divergence theorem 

0= f uAu= [ d\ud lU ) ~d l ud % u (4) 
Jn J n 



un l diU — I d l udiU (5) 
■ ><> Jn 

d i ud l u. (6) 



(ii) The source / can not be arbitrary. We integrate in Q equation J2J) to 
obtain a necessary condition for / 



= i n l d,u = Au= f. (7) 
Jan Jn Jn 

The following theorem says that J7J) is also a sufficient condition for the 
existence of solution. 

Theorem 1. A solution u to the Neumann © problem exists if and only if 
f satisfies 

// = 0. (8) 

J n 

Two different solutions differ by a constant. 

The fact that the solution is not unique in the Neumann problem does not 
affect the physics of the model that is described by these equations. Take, for 
example, Electrostatics. The electric field E l satisfies 

Ei = lU , d t E l = /, (9) 

where u is the electric potential and / the charge. If we prescribe E % ni at 
the boundary we get a Neumann boundary problem for the potential u. The 
electric field E l is invariant under the transformation u — > u + c, where c is 
a constant. We will see in section [3] that something similar happens for the 
constraint equations in General Relativity. 

We have seen that the Neumann problem has not a unique solution. If we 
include lower order terms in the operator, the Dirichlet problem will not have a 
unique solution either. For example, for some constants A > (the eigenvalues) 
the following equations have a non-trivial solutions (eigenfunctions) 

Au + Au = 0, on Q u = on dfl. (10) 



One of the main ideas in the theory of partial differential equations is that 
many relevant properties of the equations depends only on the principal part, 
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that is on the terms with highest derivatives. The previous examples show 
that uniqueness does not depend only on the principal part. Motivated by the 
Neumann problem, we write the following the two main properties of elliptic 
equations 

(i) The solutions space of the homogeneous problem (i.e., when we set the 
source / and the boundary values g equal to zero) is finite dimensional. 

(ii) The solution will exist if and only if the sources satisfy a finite number of 
conditions. 

We will see in the next sections that, under appropriate assumptions, (i)-(ii) 
depend only on the principal part of the equation and boundary conditions. 

One example of a boundary condition that does not satisfy (i) is the fol- 
lowing. 

Example 1. Let Q be the unit ball in M 3 centered at the origin. An explicit 
calculation shows that the space of solutions of the homogeneous problem 

Au = on S7, d 3 u = on dfl, (11) 

is infinite dimensional (see |27| . Chapter 1, for details). Note that the vector 
d 3 is tangential to the boundary at the points £3 = 0, x\ + x\ = 1. 



2 Second order elliptic equations 

Consider the following, second order, differential operator 

Lu = d, (a l: >(x)d 3 u + b l {x)u) + c l (x)d iU + d(x)u, (12) 

where we will assume that the coefficients are smooth functions on M. n and 
i,j = 1, • • -n. We have written the operator l)12|l in divergence form because 
it will be more suitable for the following calculations; since the coefficient a M 
and b % are smooth, this is equivalent to the standard formula 

Lu = a u (x)didju + b 3 (x)djU + d(x)u. (13) 

where b> = dja 1 ^ + V + c) and d = dib 1 + d. 

The principal part of the operator is given by the terms which contains 
only second derivatives 

l(x,d) = a ij (x)didj. (14) 

To define the symbol of L we replace in the principal part each derivative by 
the component of an arbitrary constant vector in M. n 

l(x,Z)=a i l(x)t i Z J , £eIR n . (15) 

The symbol I of L is a polynomial of order 2 in the components of £. 
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We make now the crucial assumption on the symbol. We say that the 
operator L is elliptic in J? if 

i(x, O^o y x G n, £ e K n , ^ 0. (16) 

The next important concept is the formal adjoint of L. The formal adjoint 
L f is defined by the relation 

/ vLu= / uL l v (17) 

for all u, v of compact support in J?. In this particular case we have 

L*v = dj (a ij (x)diV - c?(x)v) - b i (x)d i v + d(x)v (18) 

Note that A = A*. 

We have already seen in the case of the Laplacian that the solutions of 
the homogeneous problem play an important role; in general L and L are 
different operator and then we have two natural null spaces defined as 

Af(L) = {u : Lu = on f2 and u = on .9^} (19) 
M{L*) = {u : L l u = on Q and u = on dH}. (20) 

We can now formulate an existence result for the Dirichlct problem which will 
essentially ensures that properties (i)-(ii) are satisfied. 

Theorem 2. (i) Precisely one of the following statements holds: 

a) For each f there exist a unique solution of the boundary value problem 

Lu — f on fi, u = on dfi, (21) 

or else 

b) M{L) is non-trivial. 

ii) Furthermore, should assertion b) hold, the dimension ofJ\f(L) is finite and 
equals the dimension of J\f{L l ). 

Hi) Finally, the boundary-value problem (|21f) has a solution if and only if 

I fv = for all v G W(£*). (22) 

We will consider now the analog of the Neumann problem for L. If in the 
integration by parts given by (|17l) we allow functions u and v which are not 
of compact support, we have to include the boundary terms; and we obtain 
the following relation which is called the Green formula for the operator L 

I vL{u) - uL\v) = I vB{u) - uB^v), (23) 
J n Jon 

where the differential boundary operators are given by 
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B(u) = nja^diU + tfriiU, = n^diV - dn^. (24) 

We want to solve the following problem 

Lu = f on n, B{u) = on dfl. (25) 

As in the case of the Dirichlet problem, we define the null spaces 

Af(L, B) = {u : Lu = on 12 and S(u) = on dQ} (26) 

A/"(L*,B*) = {u : L*u = on J? and = on dJ?}. (27) 



We have the following existence result, which looks exactly the same as the 
previous theorem if we replace the Dirichlet condition by the new boundary 
condition. 

Theorem 3. (i) Precisely one of the following statements holds: 

a) For each f there exist a unique solution of the boundary value problem 

Lu = f onQ, B(u) = on dQ, (28) 

or else 

b) M{L,B) is non-trivial. 

ii) Furthermore, should assertion b) hold, the dimension of Af(L, B) is finite 
and equals the dimension of N{L i ,B t ). 

in) Finally, the boundary-value problem l|28 |l -(|24 |l has a solution if and only 
if 

[ fv = for all v £ Af(L\ B*). (29) 
J n 

We have written the boundary conditions in the form l|24(l in order to em- 
phasize that they come naturally from the integration by parts. It is possible 
to write them in a perhaps more familiar form. Define the vector /3* by 

/ 9 i = n i o ij '. (30) 

By the elliptic condition (|16|) we have /3*rii ^ 0, that is (3 l it is never tangential 
to the boundary (this excludes example [Q. I n the operator L only enters the 
symmetric part of the matrix a 4 -? , however, we have not assumed that this 
matrix is symmetric in the previous theorem. If we decompose a u = a l J + W 
where a l J — as and b^ = b^ is an arbitrary anti symmetric matrix, then 

= nj ai j +t\ T l = nj b l \ (31) 

where t % is an arbitrary tangential vector. Choosing appropriated b % and c l 
such that they do not change the operator L, we get that the function a = 
b l ni is also arbitrary. We conclude that the boundary condition B(u) = is 
equivalent to 

B(u) = (fdiU + au = 0, (32) 
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where a is an arbitrary function and [3 l is an arbitrary non tangential vector 
field on the boundary. 

Let us compare theorem [21 and |3] with the analog cases for the Laplace 
equation. We have now two operators L and L l which have two different null 
spaces (in the case when b l + c % = we have L = L* and B = B l , and then 
only one null space). There are no statements about uniqueness or about the 
elements and dimension of the null spaces. We have already seen that these 
properties depend on the lower order terms. For the particular case of second 
order elliptic operators, there exist an important tool that can give uniqueness 
and a characterization of the null space for certain kind of lower order terms: 
the maximum principle. There exist many useful versions of the maximum 
principle (see for example |14|). here we mention a particular simple case, 
which can be generalize to other situations as we will see. 

We can write the Green formula (12.'j[) in terms of a first order bilinear form 

B 

B(u,v)=f vB(u)- I vL{u)=i uB\v) - [ uL\v) (33) 
Jan Jn Jan Jn 

where 

B(u, v)= f {a lJ d u + b t u)d l v - (c l d lU + du)v. (34) 
Jn 

From this equation we deduce that u S J\f(L,B) if and only if B(u, v) = 
for all v. (One if is trivial, to see the other one, take test functions v which 
vanishes at the boundary and are arbitrary at the interior). If we assume 
b l = c % — and d < 0, then B is symmetric (i.e. B(u,w) = B(w,u)) and 
positive 

B(u,u)>0, for all u. (35) 

Moreover, B(u, u) = if and only if it is a constant and u = if d is not 
identically zero. In this case we are in a similar situation as in the Neumann 
problem for the Laplace equation: the only elements of the null space are 
the constants. More general version of the maximum principle can be used to 
prove the followings refinements of theorems [3 and [3] 

Theorem 4. Assume d < 0. Then the Dirichlet problem 

Lu = f on f2, u — g on df2 1 (36) 
has a unique solution for every f and g. 

Theorem 5. Assume d < 0, a > and not both identically zero. Let j5 % a 
vector field such that j3 l ni > on dfl. Then the oblique derivative problem 

Lu — f on Q, (3*diU + au = g on df2, (37) 

has a unique solution for every f and g. 
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In both theorems, the maximum principle can be used also to prove that the 
solution is positive if the sources, boundary values (and a in theorem (J5J) are 
positive. 

Note that in theorems |21 and [3] the null spaces for the operator and the 
adjoint have the same dimension, we will see in the next section that this will 
not be the case for more general operators and boundary conditions. 

We conclude this section with some examples. 

Example 2. The most important second order elliptic operator is the Laplacian 
on a Riemannian manifold. It is given by 

Lu = A h u = h ij DiDjU, (38) 

where h is a Riemannian metric (a lJ — h lJ ) and D its corresponding covariant 
derivative. One important example of lower order term is given by the confor- 
mal Laplacian which appears naturally in the Einstein constraint equations 

R 

Lu = A h u - —u, (39) 

o 

where R is the Ricci scalar of h a b- 

For a Riemannian metric, the principal part of the boundary condition 
B(u) has a geometric interpretation 

B{u) = n l D lU , (40) 

where we use the standard convention n l — h lJ Tij. That is, the vector n is 
now the unit normal vector with respect to the metric hij . This is sometimes 
denoted as conormal boundary condition. 

An example of lower order boundary terms is the following 

B(u) = n l D,u + Hu, (41) 

where H is the mean curvature of the boundary Q with respect to the metric 
hij. This boundary condition appears in connection to black holes (see 
and 0). 



3 Elliptic Systems 

3.1 Definition of ellipticity 

We saw in the previous section that ellipticity is a positivity condition on 
the symbol of the equation. In order to generalize this concept for systems 
of equations (this includes as particular case higher order equations) we need 
to define the symbol of a system. We can use the same idea as before, and 
define the principal part as the collection of terms which have the highest 
order derivatives. That is, consider the following differential operator in W 1 
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L(u) = a a (x)d a u, (42) 

\a\<2m 

where a is a multi-index, and the coefficients a a are N x N matrices. The 
principal part is defined as 

l{x,d)= ^ (43) 

|a|=2m 

and the symbol 

l(x,0= E a «(x)C- (44) 

\a\=2m 

The operator is elliptic if detl(x,£) ^ for every x 6 Q and £ ^ 0. This 
is the definition that appears in most text books, we will call it classical 
ellipticity (there is no general agreement on the nomenclature, in most places 
these systems are called just elliptic). This definition excludes many important 
examples, the most remarkable is perhaps the Laplace equation as a first order 
system (example^ . In order to include these cases, we need to be more flexible 
in our definition of the principal part, in particular it is important to allow 
terms of different orders in it. This is particular feature of systems which does 
not appear in higher order equations. 

It will be convenient to use a more explicit notation as the one given in 
1431) . Let u\, ■ ■ ■ , ujv be functions which depend on the coordinates x±, ■ ■ ■ ,x n . 
The operator 142(1 can be written as follows 

L ia ,(x,d)u v (x)=f ll (x), u,n= l,--- ,7V; (45) 

where L^ u are polynomials in (d\ • • • , d n ) with coefficients depending on x. 
[L^jy] is a N x N matrix, not necessarily symmetric. Note that N (dimension 
of the vectors v?) and n (dimension of R") are in general different numbers. 
Let Si, ■ • ■ , sat, ti, • • ■ , tjv be integers (some may be negative) such that 

deg(L M „) < s„+t v . (46) 

where deg means the degree of the polynomial L^ v in the derivatives d. The 
integers s M are attached to the equations and the t u to the unknowns. 

We define the principal part lfj, u (x, d) as the terms in which are exactly 
of order s^ + i„. The symbol l^{x,^) is obtained replacing in the principal 
part the derivatives by a vector £. We define the following polynomial in £ 

Z(z,0 = det(^(z,0). (47) 
The degree m of the systems is given by 

m=±deg(l(x,0), (48) 

where deg means degree in £. 

The following general definition of ellipticity was introduced in 9 
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Definition 1 (Douglis-Nirenberg Ellipticity). The system l|45|l is elliptic 
if there exist integer weights s^ and t v which satisfy 1)46( 1 and such l(x, £) 7^ 
for all real £ G MP, £ ^ 0, x S fl; where l(x, £) is given by <|47[) . 

For n = 2 we assume in addition 

Definition 2 (Supplementary condition). l(x,t;) is of even degree 2m. 
For every pair of linearly independent real vectors £ and the polynomial 
1(x,^ + t^') in the complex variable r has exactly m roots with positive imag- 
inary parts. 

Every elliptic systems in dimension n > 3 satisfies the supplementary 
condition (see PP). This no longer true for n = 2, as examplc[5]shows. A system 
that is elliptic in the sense of Dcfinitionnand satisfies also the supplementary 
condition (Definition |2J will be called properly elliptic. 

Note that the definition depends on the weights s M and t v which are not 
unique, a system can be elliptic for many different choices of weights. Also 
note that the number 2m is not related in general with the degree of the 
highest derivatives, for example for a second order system with N = 3 we 
have m — 3 (example |SJ) . The degree m is important because it gives the 
number of boundary conditions we have to impose in order to get a well 
defined elliptic problem, as we will see in the next section. 

There exists an important class of elliptic operators for which the Dirich- 
lct boundary conditions will always satisfy (i)-(ii) as we will see in the next 
section. These systems are given by the following definition. 

Definition 3 (Strong Ellipticity). The system is called strongly elliptic if 
s v = t u > and there exist a constant e > such that 



for all real £ G R n and all complex r\ 6 K . 

Note that every elliptic equation (i.e., N = 1) is strongly elliptic. Let us 
discuss some examples. 

Example 3 (Laplace equation as a first order system). This example was taken 
from 2 . Consider the Laplace equation in two dimensions 



Every equation can be written as a first order system if we introduce the 
derivatives of the unknown as new variables. That is, let u\ = d\u and u-i — 
&2U. Then we have the following system (n — 2 and N = 3) 



Re (t(^^f)>^1^5, 



■ 11 



(49) 



d\u + d\u = 0. 



(50) 



d\u x + d 2 u 2 = 0, 
d\u — u\ = 0, 
d 2 u - u 2 = 0. 



(51) 
(52) 
(53) 
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In the matrix notation 



di d 2 
dx -1 
,d 2 -1 




= 0. 



(54) 



In the classical definition, the symbol is constructed only with the terms which 
contains the highest order derivatives, in this case only with the terms with 
one derivative. Then the determinant of the symbol is 



o a a 
a o o 

6 



0, 



(55) 



and we conclude that the system is not classically elliptic. 

Take the weights t\ = 2, t% — = 1, for u, u\, u 2 and s\ = 0, s 2 — S3 
to the first, second and third equations, respectively Then, wc have 



a a 
a -1 
a -1 



a 2 + a 2 , 



(56) 



and the system is elliptic with m = 1. Another possible choice for the weights 
is the following Si —t i: with t\ = 1, t% = £3 = 0. 

Since n = 2, we have to check also that it satisfies the supplementary 
condition 

= /(e + ra) = ier + 2raa + r 2 |n 2 (57) 

then 

T± = (-cos9±isin6)\t'\- 1 \£\ (58) 

where |£| 2 = and a£i = cos ^l£ll£'l- That is, we have only one root with 
positive imaginary part. 

Example 4 (Stokes system). This example was taken from |28j . The following 
equations appear as the stationary linearized case of the Navier-Stokes equa- 
tions (see for example [33]) for the velocity u l and the pressure p of the fluid 



Au l - d l p : 

The unknowns are u l ,p, that is N = 
the matrix notation wc have 



: 0, d l u t = 0. 
4, and we will assume n 



A -d 2 
A -d 3 









u 2 






) 





0. 



(59) 
3. Then, in 



(60) 



It is clear that the system is not classically elliptic. Take t\ = t 2 = = 2, 
£4 = 1 and s\ = s 2 = S3 = 0, S4 = — 1. Then the symbol is 



fltl 2 o o 

o |ei 2 -£ 

|£| 2 -£ 

V 6 6 6 o / 
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(61) 



and we have 

Z = |£| 6 , m = 3. (62) 

Then, the system is elliptic. Another possible choice for the weights is the 
following: s$ = with £x = ^2 = *3 = 1 an d ^4 = 0. 

Example 5 (Cauchy-Riemann equation). We write the Cauchy-Riemann equa- 
tion L(u) = dgu, in terms of the real variables z — x + iy 

We have n — 2, N — 1. The symbol I = £i + i^2, satisfies 

/(£) ^ for all real £ ^ 0, (64) 

hence the system is elliptic with m — 1/2. However, it does not satisfies the 
supplementary condition because 2m = 1 is not an even number. 

Example 6. Consider the following operator in R 3 , acting on three vectors u l 

LijV? = d 3 (£u)i^ (65) 

where 

{£u)ij = 2^<9 (i Uj) + \5ijd k u k , (66) 

and /i, A are constants. Since in this case we have N = n = 3 we will use the 
same index notation for the index in the vectors it and in the coordinates of 
R 3 . 

The system (|65f) appears in elasticity (see, for example, [20]). It also ap- 
pears in General Relativity related to gauge conditions like the minimal dis- 
tortion gauge (see [33]) and in the constraint equations (see [Hi]), usually with 
the choice /i=l,A = — 2/3 which makes l|66[) trace free. 

From we deduce 

Li,** = (Ox + X)didj + fiSijA) u\ (67) 

in the matrix notation we have (A' = fi + A) 

[X'df + fiA X'd 2 d 1 X'd x d 3 
LijvP^l X'dxda X'di + fxA X'did 3 | | u 2 | . (68) 
\ X'd!d 3 X'dxd 3 X'dl+nAj 

Take Sj = ij = 1, the symbol is given by 
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I /'^ra- (69) 

and 

; = /i 2 (2 A i + A)|^| 6 , m = 3. (70) 

The operator is (classically) elliptic for /j, > 0, 2\x + A > 0. It is also strongly 
elliptic 

l^rfrf = (A + + rt%Viff > <**&W, (71) 

where e = min{/i, 2/i + A}. 

Example 7 (Einstein Constraint equations). 

There exist different ways of reducing the Einstein constraint equations to 

an elliptic systems (see, for example, the recent review U). In the standard 

approach the principal part of the system is formed with the Laplace operator 

on a Ricmannian manifold given in example |2] and the operator that has been 

discussed in example 

A particular interesting example is the one that has been recently used in 
and [7] to construct new kind of solutions. This system is not elliptic in the 

classical sense but it satisfies definition ^ f°r appropriate weights (see these 

references for details). 

Example 8 (Witten equation). The Witten equation 8aa'U A = (in the spino- 
rial notation) plays an important role in the positive mass theorem of General 
Relativity (cf. Solutions of this equation has been analyzed in |2H] and 

Em 

In the matrix notation (N — 2, and we will assume n — 3) this system is 
given by 

( d 3 di + id 2 
\di - id 2 -d 3 

The principal part, with weights ii = t 2 = 1, s± = s 2 = 0, is given by 

M^0=Q i l 3 ^/ 1 ^f 2 ) ; ^-fe 2 +£i+£i 2 ), m=l. (73) 
Then, the system is elliptic. 

3.2 Definition of elliptic boundary conditions 

For the operator L M „ defined in (|45|l we will consider boundary conditions of 
the following form 

B{x,d)i u u v = 0, l = l,---,m; u=l,---,N; (74) 

where B(x,d)i u are polynomial in d and m is given by 148fl . The order of 
the boundary operators, like those of the operators in 145|) . depends on two 




Elliptic systems 13 



systems of integer weights, in this case the system t v already attached to the 
dependent variable and a new system r\ attached to each boundary condition 
such that 

degCB,,,) <n+t„. (75) 

Note that r; can be negative and also the order of the derivatives in the 
boundary conditions can be higher than in the operator. The principal part 
bi v of the boundary operator consists of the terms in Bi v which are exactly of 
order r; +t u . 

For a given operator L, we would like to know for which boundary opera- 
tors B the solutions of the corresponding boundary value problem will satisfies 
(i)-(ii). The answer to this question is given by the following definition, as we 
will see in the next section. 

Let xo a point on dfl and let n l the outer normal to fl. We consider the 
constant coefficient problem 

l ia ,(x ,d)u v = Q, (76) 
h„(x ,d)u u = 0, (77) 

on the half plane (x l — x ) ■ rii < with boundary (x l — x l )ni = 0. 

Definition 4 (Complementing condition). We say that the complement- 
ing condition holds at xq if there are no nontrivial solutions of <|76|) - (|77|) of 
the following form: 

u u {x)=v u {rj)e i ^ xS - xi o) (78) 

where £ is a any nonzero, real, vector which satisfies S^ni = 7 v(rj) tends 
to zero exponentially as rj — > — oo and the coordinate r\ is defined by rj = 
(x j - <4>) n j- 

In the literature, these conditions are also called Lopatinski- Shapiro condi- 
tions or covering conditions (see [2] and [SHj). Let us study some examples of 
boundary conditions. 

Example 9 (Boundary conditions for the Laplace equation.). Consider solu- 
tions of the form l|78|) for the Laplace equation Au = 0. We chose coordinates 
in R" such that n — x n , n % — 8 l n . Then, all the solutions of this form are given 

by 

where £ satisfies £ n = 0. 

We consider different boundary conditions on the plane x n = 0. For the 
Dirichlct condition u{x„ = 0) = we get 

u = e^ x * = 0, (80) 

since this is not possible there exist no solution of this form which satisfies 
the Dirichlet conditions. Hence, the Dirichlct boundary condition satisfies the 
complementing condition. 
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For the Neumann condition we have d Xn u = at x n = 0, this implies 
£ = and then the solution will not decay at infinity. Hence, the Neumann 
conditions satisfied also the complementing conditions. 

Take the oblique derivative boundary condition (3 l diU = at x n = 0. This 
implies 

i(Af)=0, = 0. (81) 

If P n ^ then |£| = 0, an the complementing condition is satisfied. This was 
the case studied in section [21 On the other hand, if j3 n = (like in example 
then the complementing condition is not satisfied since we can always chose 
a vector £ such that = and we will get solutions of the form (f7H|l . 

Consider now the following interesting example studied in 15 . At x n = 
we impose the boundary conditions 

Su = 0, (82) 

where 

Su = d\u + • • • + dl_ lU , (83) 

is the Laplacian in n — 1 dimension. From (|82|) we deduce the |£[ 2 = and then 
it satisfies the complementing conditions. It is also clear that S k u = where 
fc, is an arbitrary natural number, satisfies the complementing condition. Note 
that in this cases the boundary operator has derivatives of higher order than 
the Laplace operator. On a Riemannian manifold, these conditions can be 
written in geometric form where 6 is the intrinsic Laplacian on the boundary. 
Another interesting condition which also satisfies the complementing condition 
is the following 

5u - n l d z u = 0. (84) 

In this case, integrating by parts, it is easy to show that the only solutions of 
the homogeneous problem are the constants 

= / uAu = f un l diu - \ diud l u (85) 
J (2 Jan J (2 

= (p uSu — / diud l u (86) 
JdC2 Jn 

= - I \du\ 2 - [ d lU d l u, (87) 
Jon Jf2 

where du denotes the gradient intrinsic to the boundary. 

Example 10. Consider the operator discussed in example Integrating by 
parts we get 

B(u,v) = - [ v l L ljU 3 + <£ (Su) ij n i v j (88) 
Jn Jan 

where 

B(u,v) = [ (SuY^iVj. (89) 
Jn 
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We can write the integrand in B(u, v) in the following form 

{£ uyi&vj = |(£u)y(£w)« + (A + ^)d k u k d lV l , (90) 

where {Cu)ij is the trace free part of d^Uj), that is 

2 

(£u)ij = 2d(iUj) - -S i:j d k u k . (91) 

Note that B is symmetric B(u, v) = B(i>, u). Using this and equation 188(1 we 
get the following Green formula 

/ v t L ij u j — tfLijV 3 = (h (£w) i:) -nV - (£v)ijTi l vP . (92) 
Jn Jan 

This is analogous to the Green formula for second order equations (J2HJ- For 
simplicity we have not included terms in non divergence form in the operator, 
that is why we have L — L l and B = B f in H92fl . these extra terms can be 
handle in the same way as in section [21 

The boundary integral in the Green formula Ij92(l suggests that two natural 
boundary conditions are the Dirichlct 

u l = on <9J7, (93) 

and the analog to the Neumann boundary condition 

(£u) lj n j = on dtt. (94) 

We want to prove that these boundary conditions satisfy the complementing 
conditions. We will assume that /i > and 2fi + A > 0, that is, the operator is 
elliptic as we have seen in exampleEl We will make also an extra assumption: 
3A + 2/i > 0; this implies that the integrand H90(l is positive. Moreover, if 
3A + 2^i > 

B(u, u) = d {l u j} = 0, (95) 

that is u is a Killing vector. If 3A + 2/i = 0, then 

B(u,u) = ^f=> (Cu) i:j = 0, (96) 

then u is a conformal Killing vector. In flat space, we know explicitly all the 
Killing and conformal Killing vectors. Hence, we have a characterization of the 
null spaces for these boundary conditions. The Killing and conformal Killing 
are the analog of the constants for the Neumann problem for the Laplace 
equation. 

Assume we have a solution u of the form 1)78(1 . Chose Cartesian coordinates 
such that 77 — X3. Let L\ = 27r/£i and L2 = 2tt/^2- Take as domain the 
infinite cubic region X3 > 0, < x\ < L\, < X2 < -^2- For this domain 
we use equation (18811 for u = v. We want to prove that, on this domain, the 
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boundary integral in (|88|l vanished if we impose either l|93[l or (|94|l . Using these 
boundary conditions we get that the integrand vanishes on the face x 3 = 0. 
The integrand also vanishes on the face £3 = 00 because the solution, by 
hypothesis, decay at infinity. On the other faces the integrand does not vanish. 
However, because of the choice of L\ and L2, we have that the integrand of 
opposite faces are identical. Then, the sum of the boundary integrals vanished 
because the normal is always outwards. We conclude that B(u,it) should 
vanish. But there are no Killing or conformal Killing vectors which decay to 
zero at infinity. Hence the complementing condition is satisfied. 

Example 11 (Boundary conditions for the Stoke system). If we multiply equa- 
tions 023 by u l and integrate by parts we get 

= - J d kUl d k u l + <j> n k (u l d k u t - u kP ). (97) 

Using this equation and a similar argument as in the previous example it is 
possible to show that the boundary conditions 

u l = on dn, (98) 

and p unprescribed, satisfy the complementing condition (see for example 

ESI)- 

Example 12 (Dirichlet boundary conditions for strongly elliptic systems). 

Assume that the system is strongly elliptic, this implies Sj = tj = > 0. 
The Dirichlet boundary conditions on df2 are given by 

(n%)% = 0, 9 = 0,--- -1, j = l,---,N; (99) 

when tj — 0, Uj goes unprescribed. 

It can be proved that for every strongly elliptic system, the Dirichlet con- 
ditions (|99fl satisfy the complementing condition (see [2]). 

In the case equations (N = 1) of order 2m (ti = m) these conditions 
reduce to 

{n l di) q u = 0, q = 0,---,m-l. (100) 

In particular, for second order equation (m = 1) we have u = at the bound- 
ary. That is, we recover the familiar Dirichlet condition studied in sections ^ 
and |21 In example we have t[ = 1, then q — and the Dirichlet conditions 
is just 

u j = on df2. (101) 
As an example of a higher order equation, we have the biharmonic equation 

AAu = f, (102) 
the Dirichlet conditions are given by (N — 1, m = 2) 



u = 0, n l diU = on dQ. 



(103) 
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Example 13. In the following example (taken from |28p. the complementing 
condition is not satisfied. Consider the following problem in R 2 , where the 
boundary is the line x 2 = 

AAu = on Q Au = d 2 Au = on dQ. (104) 

For every £ £ E the function 

u(x,y) = e^ xl ~^ X2 (105) 

is a solution. 

Example 14- We have seen that for strongly elliptic system the Dirichlet 
boundary conditions satisfy the complementing condition. This is not true 
for general elliptic systems. In following example (discussed in |2S1) we show 
that there are elliptic systems for which the Dirichlet problem is not well 
defined. 

Consider the system (N = n = 2) 



3\ - d\ (ui 
29i5 2 d\ - dl) U 2 



= 0. (106) 



The symbol is given by 



l - = K^ i- ti)' l = ® + ® - (107) 

Then, the system is elliptic in the classical sense. This system can be written 
in the complex form, z = x\ + ix 2 , w = u\ + iu 2 , as 

dlw = i(o»i + id 2 fw = 0, (108) 

for which the general solution is clearly 

w = f(z) + zg(z), (109) 

where / and g are arbitrary functions of z. We observe that all solutions of 
the form 

w = f( z )(l -zz) on \z\ < 1 (110) 

with arbitrary analytic /, vanish on the boundary of the unit disk. Thus the 
problem of finding a solution with Dirichlet boundary condition is not well 
defined. 

Example 15 (Boundary conditions for the Witten equation). In the Witten 
equation studied in example [S] we have m = 1, that is, we can only impose 
one boundary condition. Consider the following boundary condition 

U! = 0, on dQ, (111) 
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and u 2 goes unprescribed. This condition has been studied in as an inner 
boundary condition for black holes in the positive mass theorem. We want to 
prove that it satisfies the complementing condition. We can explicitly calculate 
all the solutions of the form (|78|l of the equations l|72|) 

u v = e^ x *v u (x 3 ), v u = A"e^ X3 , (112) 

where A" are constants such that A 2 /Ai — (i£x + ^ 2 ) I ^ I 1 an d we chose 
coordinates such that r\ = X3, £3 = 0. There is no solution of this form that 
satisfies ui(x3 = 0) = and then the complementing condition follows. 

Example 16 (Stationary solutions of Einstein equations). 

In the presence of a timelikc symmetry, Einstein equations can be reduced 
to an elliptic system. Moreover, the inner boundary conditions satisfy the 
complementing condition. This result was proved in |3(J) and it was used to 
prove an existence result for the non linear problem. See also |2.'jj for a different 
kind of boundary conditions for the static case. 

3.3 Results 

In order to present a general result for properly elliptic systems with boundary 
conditions that satisfy the complementing condition we need to reformulate 
in a more precise way properties (i)-(ii). For a given operator L and boundary 
operator B we consider the operator A defined as A(u) = (L(u), B(u)). This 
operator will act on appropriate Hilbert spaces Hi and 7i. 2 , A : Hi — > 7i 2 . In 
analog way as we did in (1261) we define the null space M(A) as 

Af(A) = {u G n x : A(u) = 0}. (113) 

The range of A is defined by 

K(A) = {weH 2 :3u€H u A(u)=w}, (114) 

and the complement of the range is given by 

K^{A) ={weH 2 : (Au,w) H2 = for all u G Hi}. (115) 

where (•, •) denotes the Hilbert scalar product. 
We can write now properties (i)-(ii) as follows 

(i) M(A) has finite dimension. 

(ii) 72. (A) has finite dimension. 

An operator which satisfies (i)-(ii) is called a Fredholm operator, (we have 
assumed that A is bounded, otherwise in (ii) we need to impose that 71(A) is 
closed, see ^Sl an d [TS|1. 

We have the following general result (we only sketch the statement, for 
details and proofs see ^Tj^CHI an d also EE]) 
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Theorem 6. If the system L is properly elliptic in D and the boundary con- 
ditions satisfy the complementing condition for every point of dfl, then the 
operator A(u) = (L(u), B(u)) is Fredholm. 

We have seen that the dimension of M (and hence uniqueness) is not 
invariant if we add lower order terms to the operator. One of the consequence 
of theoremElis the existence of an invariant for elliptic problems: the Fredholm 
index. This number is defined as 

I = dimTV(A) - dim ft" 1 (A). (116) 

It can be proved that the index I is stable under perturbation, in particular 
it does not depend on the lower order terms. 

In section|5]we have used the Green formula to construct the formal adjoint 
operator L l and its corresponding boundary operator B . In this case we can 
define A\u) = (L*(u), and it can be proved 1 that N(A f ) = K X (A). 

That is, the boundary value problems considered in theorem [21 and [21 have 
7 = 0. In fact these theorems also show that the index does not depend on 
the lower order terms in this particular case. 

Boundary conditions which come from a Green formula are called normal 
boundary conditions. The advantage of them is that we have a characterization 
of TZ ± (A) through the formal adjoint problem, and then we can in principle 
compute the conditions that the sources should satisfy in order to have a solu- 
tion. General results for normal boundary conditions for higher order elliptic 
equations can be found in |19| |32j . Since these boundary conditions come 
from an integration by parts, the order the boundary operators will be al- 
ways less than the operator itself. We have seen that this is not necessary the 
case for general elliptic boundary conditions that satisfy the complementing 
condition. For the general case, we will not have a characterization of TZ ± (A). 

We have seen that the Dirichlet boundary conditions satisfy the comple- 
menting condition for strongly elliptic systems. Using this fact, general exis- 
tence results for the Dirichlet problem can be proved (see (22)- Moreover, it 
can be shown that the index is always zero in this case. 

Finally, we want to present an existence result for the operator considered 
in example [5] that can be deduced from the general theorem[Hl(see |3S|) In this 
case, we have a Green formula and then we have normal boundary conditions. 
The following two theorems are the analogous of theorem [21 and [31 

Theorem 7. Let Lij be given by (1651) with [i > 0, 2/i + A > 0, 2[i + 3A > 0. 

Then, for every smooth f J and g % , there exist a unique, smooth, solution u l 
of the Dirichlet problem 

L l3 u l = f on Q, u l = g l on dQ. (117) 

1 It is important to note that for any bounded (or unbounded with dense range) 
operator A we can define the Hilbert adjoint A' . This is not related, in general, 
with the formal adjoint A . However, when we have a Green formula, it is possible 
to prove that in fact A 1 — A' (see theorem 8.4 of |19|. |S] and also |35| 
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We have seen that all solution of the homogeneous problem satisfy 
0. that is v is a Killing or a conformal Killing vector. Uniqueness in this 
theorem follows because there exists no Killing or conformal Killing vector 
which vanishes at the boundary. 

Theorem 8. Let Lij be given by (|65[) with /i > 0, 2/i + A > 0, 2/i + 3A > 0. 

Consider the boundary value problem 

L %j u l = f on Q, (£u) ljU l = on dO. (118) 
This problem has a solution if and only if 

/ fiV 1 = for all v l such that £vij — 0. (H9) 

If Ui andu2 are two different solutions, then the difference v = Ui — u 2 satisfies 

= o. 

In the case of the Einstein constraint equations, the previous theorems can be 
used to prove existence of solutions of the momentum constraint (see US]). In 
this case the physical quantity is the second fundamental form Kij which is 
given by 

K ij =Qi j -(£v)i j , (120) 

where Q is an (essentially) arbitrary tensor. Then, as in the case of the Neu- 
mann problem for the Laplace equation, the lack of uniqueness in theorem [S] 
will not affect K«j. 



4 Final Comments 

In order to check if a system of equations is elliptic, we should first prove that 
the principal part of the operator satisfy definitions ^a-nd [21 If the system is 
non linear, we should consider the corresponding linearized problem. Then we 
should prove that the boundary operators satisfy the complementing condition 
(definition^. This can be complicated. There exist equivalent formulations 
of this condition (see for example |38j ) . some of them can be more suitable 
for specific problems. It is also important to know if the boundary conditions 
come from a Green formula (normal boundary conditions). The Green formula 
can be used to prove that the complementing condition holds (as we have seen 
in the examples). Moreover, for the case of higher order equations there exist 
general results that can be used (see 

We have only discussed linear elliptic systems. In the non linear case, there 
are no general existence results like theorem [H] For non linear second order 
equations a good reference is |14j and for non linear systems |13| and |12j . A 
related issue that was not discussed here is regularity. We have assumed that 
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the functions and the boundary are smooth. Regularity properties are crucial 
for non linear systems, see H3 and and references there. 

In sectional we have followed JJ. Good references for this section are also 
|10|. [Tl| and [52]. For sectional an introductory book is [2E]; more advanced 
material can be found in and for a complete discussion see ^B], \To\ 
and [3HJ 
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